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MIDCOURSE GUIDANCE FOR RETURN FROM THE MOON TO A 


GEOGRAPHICALLY FIXED LANDING SITE 

By John D. McLean and Luigi S. Cicolani 
Ames Research Center 

SUMMARY 

This report describes a midcourse guidance system for return from the 
moon to a safe landing at a specific geographical site. The design is based 
on linear perturbations about a reference trajectory. The vehicle's ability 
to maneuver within the atmosphere is used to reduce the midcourse corrective 
velocity requirements. 

The system is compared statistically on the basis of a digital computer 
simulation with one that uses fixed^time-of -arrival method of guidance. The 
trajectory estimation method used for both systems is that described in NASA 
TR R-135* The transearth injection errors^ observation errors^ velocity cor- 
rection mechanization errors^ and velocity correction measurement errors are 
specified by statistical distributions considered to be realistic in terms of 
present day capabilities. 

The fixed-landing-site guidance system requires substantially less cor- 
rective velocity than a fixed-time -of -arrival system, but it is more sensitive 
to errors in the final velocity correction than is the fixed -time -of -arrival 
system. This sensitivity can be effectively compensated for by proper schedul- 
ing of observations and velocity corrections or by the use of a rocket engine 
having different error characteristics for midcourse corrections. 


INTRODUCTION 


The objectives of midcourse guidance for the return phase of a lunar 
mission are defined for this study as: 

1. To place the vacuum perigee of the return trajectory at the center of 
the entry corridor. 

2. To land the vehicle at a specified geographical landing site. 

3 . To reach this landing site without crossrange maneuvering during 
entry. 

The third objective arises because the crossrange maneuvering capability for 
Apollo -class entry vehicles is much more limited than the downrange maneuver 
capability (ref. l) . 

It is recognized that because of errors, from sources to be discussed 
later, the above objectives cannot be met exactly. However the 



fixed-time- of -arrival guidance system used in earlier studies at Ames (ref. 2) 
and elsewhere does not explicitly impose these objectives. (This system 
guides the vehicle to the position on a precomputed reference trajectory cor- 
responding to the time of reference perigee.) It is true that the fixed-time- 
of-arrival system satisfies the desired constraints within acceptable 
tolerances for the range of initial errors studied in reference 2. Vacuum 
perigee differs from the reference value (center of the entry corridor) by 
small amounts and the desired landing site can be obtained with modest cross- 
range and downrange maneuvering during the entry flight. On the other hand^ 
there may be some advantage in a guidance system which recognizes these 
c ons t r aint s expli c it ly . 

This report presents the development of a method for guidance to a fixed 
landing site. The perfoimance of this system for a sample return trajectory 
is then compared statistically with that of a f ixed-time-of-arrival system on 
the basis of a digital simulation. 


NOTATION 


C covariance matrix of velocity corrections actually made 
D declination 

G 3x6 matrix from guidance equation 

H matrix of partial derivatives of observed angles with respect to 
Cartesian coordinates 

I unit matrix 

kg multiplying factor for initial covariance matrices 

M 3x3 matrix from guidance equation 

P covariance matrix of errors in trajectory estimation 

Q covariance matrix of observation errors 

R magnitude of R 

R position vector 

r vector of small deviations from reference position 

rms root mean square 

RA right ascension 

S covariance matrix of errors in making velocity corrections 
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V 

V 


V 

W 


X 


a 

7 

6 yO 

9 

0 

U) 


covariance matrix of errors in measuring velocity corrections 
time 

unit vector 
magnitude of V 
velocity vector 

vector of small deviations from reference velocity 
corrective velocity increments 

covariance matrix of deviations between actual and reference 
trajectories 

state vector (6x1 matrix of vehicle’s position and velocity 
deviations from reference) 

azimuth angle from north of orbital plane at landing site 

one-half the subtended angle of earth or moon 

pointing errors in corrective velocity 

4x1 matrix of deviations in arrival parameters 

entry range angle 

transition matrix 

matrix of partial derivatives of Cartesian position and velocity with 
respect to arrival parameters 

earth’s angular velocity 


Superscripts 

T transpose of a matrix 

( ) 3x1 matrix (or vector) 

(*) derivative with respect to time 

Subscripts 

a actual 

d desired 
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E atmospheric entry 

i^m integers 

L landing site 

p perigee 

rms root-mean-square value 


ANALYSIS 


The development of the method of guidance to a fixed landing site follo-ws 
the approach used in reference 3 for finding reference trajectories vhich 
return from the moon to a specified geographic site. Such a reference tra- 
jectory distorted for ease of illustration is pictured in sketch (a) . The 



pointy represents either the transearth injection point or any point along 
the trajectory. The trajectory is divided into tvo phases; orbital and entry^ 
entry being defined as a radius of 65 OO km (about 400^000 ft altitude). 

The point of view is taken that the trajectory travels backward in time 
from the landing site to point P. Hence ^ the orbital portion can be speci- 
fied by the vacuum perigee position vector Rp, perigee velocity vector Vp^ 
and time of perigee^ tp. If the trajectory is to terminate at the landing 
site^ these quantities must satisfy certain functional relationships. These 
are; 

Rp = Rp 

Vp = \ (9,cc,tL,DL) 

^p = ■*^L “ 0-009339 - 0.00254 


( 1 ) 

( 2 ) 


k 



where cp is the angle from entry to landing; the. azimuth angle, a, is mea- 
sured from the meridian through the landing site to the orbital track; is 

the declination of the landing site, and is the time of landing. Equa- 

tions (l) are geometrical relationships (derived in appendix A) which implic- 
itly restrict the crossrange maneuver to zero. Equation (2) contains the 
dynamical relationship between the entry angle, cp, in radians and the time in 
mean solar days required to traverse it. This equation is of empirical origin 
(see ref. 3 ) • Thus the complete trajectory may be specified in terms of the 
"arrival parameters" cp, a, t^, Vp, Rp^ and 

Assume that a reference trajectory from the moon to the landing site has 
been defined. The point P in sketch (a) will now be regarded as lying on 
the reference trajectory somewhere in cislunar space. Because of injection 
errors, the vehicle arrives at point P’ at time t instead of at P (see 

sketch (b)). Likewise, the 
velocity will generally be 
such that unless a velocity 
correction is made before 
entry the vehicle will not 
land at the desired site. 

If a velocity correction 
is to be made at time t, 
then the function of the 
fixed-landing-site guidance 
system is to compute a new 
trajectory which leaves P’ 
at time t and lands at the 
landing site at time t^. 
Because of the earth's rota- 

Correcfed orbit for fixed lending site tion during the interval 

(t{i - tjj) the landing site 
will have moved to the point 
L' . The corrective velocity increment will be the difference between the 
vehicle's actual velocity and that of the corrected orbit. 

For comparison the corrected trajectory which would be produced by the 
fixed-time-of-arrival method is shown as a dashed line in sketch (b). It 
intercepts the reference at P^, corresponding to perigee on the reference 
orbit. Note that if the reference perigee point lies outside the orbital 
plane, a plane change will be required. On the other hand, if Rp and 
are the only arrival parameters constrained, the fixed- landing- site guidance 
system will require a plane change only if the inclination of the incoming 
orbit is less than (see ref. 3 ) • 

Finally, it should be pointed out that the basic methods used in 
developing the fixed- landing-site guidance system can be used for landing on 
any rotating celestial body, provided the "entry" phase of the trajectory can 
be suitably approximated. This statement is true even if the destination body 
has no atmosphere (here "entry" would correspond to powered descent), but it 
would probably be pointless to apply the method when the rotation rate of the 
body is extremely low. 
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Fixed-Landing-Site Guidance Equation 


For the computation of a new trajectory from P' to the landing site^ 

% and Dl are constrained to the reference values . This leaves four arrival 
parameters (cp^ cl, t^^ and Vp) which may be changed from the reference values 
so that the trajectory will originate at P' instead of P. Since only three 
quantities are required to specify the position of P'^ there is an infinite 
family of possible corrected orbits. A unique trajectory may be selected from 
this family either by constraining an additional arrival parameter or by sat- 
isfying some other condition. For this study we require magnitude of the 
corrective velocity to be a minimum. 

The procedure for finding the reference trajectory is a complex iterative 
one^ so the problem of finding the corrected trajectory has been linearized in 
terms of small perturbations around the reference trajectory. The resulting 
linearized guidance equation is derived in appendix A and may be written as 

= (I - uu'^)(Mr - v) ( 3 ) 

where v^ is the desired corrective velocity vector^ u is a unit vector^ M 
is a 3 X 3 matrix and r and v are vectors of position and velocity 
deviation from the reference trajectory. 

Comparison with fixed time o_f_ ar rival.- The f ixed-time-of-arrival 
guidance equation in reference 2 is given by 

^d ~ <l>2“^ir-v (4) 


where O2 and are submatrices of the state transition matrix which relates 
position and velocity deviations at the time of reference perigee to those at 
the time of the correction. 

If the vector (Mr - v) in equation (3) is defined to be the nonoptimiun 
correction^ then both the f ixed-time-of-arrival and the nonoptimum fixed- 
landing-site guidance corrections can be expressed in the form (Ar - v) . It 
can be seen that the magnitudes of both corrections are dependent on the 
direction of Y , but the correction will be zero only when Ar is equal to v. 

Now consider the multiplying factor (l - uu"^) in the landing-site 
guidance law. Since 


(l - = 0 


any component of the nonoptimum correction lying in the u direction will be 
eliminated by this multiplying factor. Because of this fact^ it is to be 
expected that the magnitude of the optimum correction for fixed- landing-site 
guidance will be more sensitive to the directions of r and v than is the 
case for f ixed-time-of-arrival guidance. Thus^ from a statistical point of 
view^ the mean-square velocity correction for fixed- landing-site guidance will 
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be more sensitive to the directional properties of the distribution of 
deviations from the reference trajectoiy. 


DIGITAL SIMULATION 


The fixed- landing-site guidance equation was compared on a statistical 
basis with that for fixed time of arrival by means of a digital computer 
simulation of the two systems under identical conditions . 

Four sources of error were considered: 

(1) Transearth injection errors 

( 2 ) Errors in knowledge of the vehicle position and velocity (trajectory 
estimation errors) 

(3) Errors in mechanization of the desired velocity correction 

(4) Errors in measurement of the velocity correction actually made 

The injection and estimation errors were assumed to be identical for each sys- 
tem at injection time and to be described statistically by a knovn covariance 
matrix. The particular covariance matrices used in the study will be 
discussed in the section on results. 

During the course of the flighty celestial observations are made in order 
to reduce the trajectory estimation errors. The resulting data are processed 
by an optimal filter (see ref. 4) which provides the best estimate of the 
spacecraft's position and velocity together with the covariance matrix of 
errors in that estimate. The type of celestial observations and the measure- 
ment error model assumed are identical with that used for the standard case in 
reference 2 , as is the error model for item (4) above. These error models are 
described briefly in appendix B. While it is necessary to separate the errors 
of items (2) and (4) in the mathematical analysis^ they will be lumped 
together under "trajectory estimation errors" in the remainder of this 
discussion. 

The error model for the velocity correction^ the one used in reference 5> 
is described in detail in appendix B. It is shown in the appendix that for 
small velocity corrections (of the order 1.0 m/sec) the vector error in the 
correction is due almost entirely to rocket engine cutoff error. This error 
lies nearly in the direction of the commanded velocity correction and has an 
rms value of 0.1 m/sec. On the other hand^ as the magnitude of the commanded 
correction increases^ the error normal to the commanded direction increases 
until the components of the error along and normal to the commanded correction 
are nearly equal. (The rms value of both components in this case depends on 
the magnitude of the correction.) 
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Schedule of Ohservations and Velocity Corrections 


The schedule of observations and velocity corrections used is considered, 
on the basis of other Ames studies, to be nearly optimum for the f ixed-time-of- 
arrival system. This schedule is outlined briefly as follows: 

(1) Starting a half hour after transearth injection, ten observations 
are made followed by the first velocity correction, all at half-hour intervals. 

(2) Starting 24 hours after transearth injection, nine observations are 
made followed by the second velocity correction, all at one-hour intervals, 

(3) Starting 34 hours after transearth injection, 15 observations are 
made followed by the third and final velocity correction, all at one-hour 
intervals . 

This last correction occurs about 8.6 hours before entry. 


Computation of Statistical Information 

The purpose of the digital computer simulation was to determiine how well, 
and at what cost in corrective velocity, each guidance system fulfills the 
requirement of returning the spacecraft from the moon on a trajectory from 
which a satisfactory landing can be made. The two guidance systems are then 
compared on the basis of: (l) the rms values of the velocity corrections at 

the three different correction times and of the total correction, and (2) the 
rms deviations from the reference values of the "landing parameters," perigee 
altitude, entry range angle, crossrange adjustment, and time of landing. 

During the orbital phase of the trajectory, the statistics of the 
velocity corrections and deviations from the reference trajectory were com- 
puted using the linear methods described in appendix C. However, it was nec- 
essary to use the Monte Carlo method for the entry phase of the flight. The 
computations involved for the entry phase are outlined in appendix D. 

The first velocity correction attempts, on the basis of the estimated 
state vector at the time of the correction, to eliminate the effects of the 
injection errors. It is only because the estimate of the vehicle's position 
and velocity and the mechanization of the velocity correction are imperfect, 
and hence sources of error, that the second and third corrections are needed. 

It is shown in appendix C how the portions of the second and third corrections, 
for these sources of error, can be separated. 

Finally, in the case of fixed- landing-site guidance the deviations from 
the reference values of entry range angle and time of landing arise from two 
sources and can be separated accordingly (see appendix C) . Some of these 
deviations, as in the case of fixed time of arrival, result from errors in 
trajectory estimation and velocity correction mechanization and will be 
referred to as errors . The remaining deviations result from changes required 
by the guidance law and will be referred to as entry range and landing time 
adjustments , 
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RESULTS AND DISCUSSION 


Data from the eight cases outlined below provide a comparison between the 
fixed- landing-site and f ixed-time-of-arrival systems with regard to corrective 
velocity requirements^ terminal accuracies and the effects of the three major 
error sources (injection errors^ observation errors^ and velocity correction 
mechanization errors) . These data are intended to indicate the areas which 
must be considered for a particular application of the guidance law. 

Two different covariance matrices of transearth injection errors were 
used. The first is considered^ on the basis of unpublished work^ to be realis- 
tic for an Apollo-type mission. This distribution has mean -square values of 
about 0.6 kra^ and 3*5 m^/sec^ in position and velocity^ respectively. However^ 
the errors are cross correlated and the complete statistical description is 
given by the covariance matrix presented in table 1 . This distribution was 
used for most of the data and will be referred to as the ^standard" 
distribution. 

The other covariance matrix of injection errors is of interest mainly 
from a theoretical point of view. This matrix is diagonal for Cartesian coor- 
dinates^ with mean -square values of 3 km^ and 3 m^/sec^ in position and veloc- 
ity. The mean-square errors are equal for each Cartesian direction^ so the 
distribution will be referred to as spherical. 

The two covariance matrices of injection errors were multiplied by a 
scalar^ kg, in order to assess the influence of the magnitude (as opposed to 
direction) of the initial error distribution. It is shown in reference 2 
and is confirmed here that the effects of the magnitude of the injection 
errors on the estimation errors are of minor significance. 

The effects of velocity correction mechanization errors were evaluated by 
comparison of data resulting from the use of the standard error model des- 
cribed in appendix B with the data resulting when the error was assumed to be 
zero . 


With this background the eight cases for which data are to be presented 
can be summarized as follows: 


Case K 


la 

1 

2a 

10 

3a 

100 

lb 

1 

2b 

10 

3b 

100 

Ic 

1 

3c 

100 


Velocity 

correction 

errors 

None 

None 

None 

Standard 

Standard 

Standard 

Standard 

Standard 


Injection 

errors 


Standard 

Standard 

Standard 

Standard 

Standard 

Standard 

Spherical 

Spherical 
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Corrective Velocity Requirements 


The ratio of the rms error in estimating the vehicle's position for case 
3 a to that for la is plotted in figure 1 as a function of the number of obser- 
vations. Note that just before the first velocity correction the estimation 
error for case 3 n is only about I/3 larger than that for case la even though 
the initial rms errors differed by an order of magnitude. This increase in 
estimation error vith initial error becomes progressively smaller as more 
observations are made^ and the results are substantially the same for velocity 
estimation errors. Thus^ it can be seen that the errors remaining to be cor- 
rected after the first velocity correction are mainly the result of observation 
errors andj, where applicable^ of velocity correction mechanization errors. The 
rms values of the individual velocity corrections and the total corrective 
velocity are presented in table II. Likewise^ the ratio of the total rms cor- 
rective velocity for fixed landing site to that for fixed time of arrival is 
presented. As expected on the basis of the data in figure 1 ^ the first veloc- 
ity correction increases linearly with^/k^ while subsequent corrections show 
only a small effect from initial errors. 

Total corrective veloc ity requirements,- The total rms corrective velocity 
for fixed- landing-site guidance ranges from 5^ percent of that for fixed time 
of arrival in case la to 6 k percent in case 3 n. Including the velocity cor- 
rection mechanization errors changes the result for the first case to 5^ 
percent but has negligible effect on the third case because of the dominant 
influence of injection errors. These data indicate that the fixed- landing- 
site guidance system requires substantially less corrective velocity than 
fixed time of arrival for a realistic set of injection errors. However^ for 
cases Ic and 3c the corrective velocity requirements for fixed- landing- site 
guidance increase to 70 and 89 percent of those for fixed time of arrival^ 
indicating that the initial distribution of errors must be considered carefully 
in comparing the two systems . 

The importance of the directional distribution of the deviations being 
corrected for is demonstrated by comparison of the data in table II for cases 
Ic and 3 c with those for cases lb and 3 b. The spherical distribution of injec- 
tion errors requires an initial correction for fixed- landing- site guidance 
which is about 95 percent of that for fixed time of arrival as compared to 
69 percent for the standard distribution. The remaining corrections are in 
the same ratio as for the standard distribution^ since the first correction 
eliminates most of the effects of initial errors. 

Penalties from veloci ty correction mechanization errors.- Table III lists 
the rms values of the portions of the second and third velocity corrections 
which result from errors in making the previous corrections. The portion of 
the second correction due to errors in the first increases with initial errors 
because the magnitude dependent portion of the velocity correction mechaniza- 
tion error in the first correction becomes significant for cases 2b and 3b. 

On the other hand^ since the second correction is small the rms error in the 
correction is essentially constant and the rms corrective velocity it requires 
at the time of the third correction is constant . 
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Fixed- landing-site guidance requires only about half as much corrective 
velocity at the time of the third correction because of errors in the second 
correction as does fixed time of arrival even though the rms mechanization 
errors are the same. However^ from the data in table II (cases la^ 2a^ and 
3a) it is seen that the ratio of the corrections caused by trajectory estima- 
tion errors at the time of the second correction was about l/5* This differ- 
ence is due largely to the fact that for small corrections^ the mechanization 
error lies approximately along the commanded velocity correction. For fixed- 
landing-site guidance this means the error has no component along the unit 
vector u while one would not expect this to be the case with the estimation 
errors, (it should be emphasized here that the magnitude of the third correc- 
tion is rather insignificant and it is used only to demonstrate a directional 
sensitivity which could be important in other applications of the fixed- 
landing-site guidance scheme.) 


Terminal Deviations 

Perigee altitude .- The rms deviations from reference perigee altitude are 
summarized in table IV. Only the errors in the final correction and in the 
observational data have an appreciable effect on deviations from the desired 
terminal conditions . Therefore the deviations for the first three cases and 
the last three ^ respectively^ are nearly constant. In addition^ the devia- 
tions for the two guidance systems are nearly equal in the absence of velocity 
correction mechanization errors . 

The data in table IV also indicate that perigee altitude is much more 
sensitive to errors in the final velocity correction for fixed-landing-site 
guidance than for fixed time of arrival. It was pointed out earlier that for 
velocity corrections of the order of 1 m/sec ^ the mechanization error results 
almost entirely from the rocket engine cutoff error. Therefore^ if the final 
velocity correction is small^the rms perigee altitude error will be determined 
almost entirely by the magnitude of the cutoff error and the time of the final 
correction. The perigee altitude error due to errors in the final correction 
is much larger for the landing-site-guidance system than for f ixed-time-of- 
arrival guidance. For this reason the use of a rocket engine with a smaller 
cutoff error would improve the accuracy of the former scheme relative to 
fixed-time -of -arrival guidance. It can be shown analytically (on a two-body 
basis) that perigee altitude becomes less sensitive to velocity changes as the 
time to go (to perigee) decreases. Therefore^ since the mechanization error 
is nearly constant^ the perigee error could be reduced^ at the expense of more 
corrective fuel^ by delaying the final correction. This approach will be 
discussed in more detail later. 

Crossrange .- The rms crossrange deviations for the two systems are 
presented in table V. (See appendix D.) The fixed-time-of- arrival system 
does not attempt to eliminate the necessity for crossrange adjustment during 
entry so that the deviations for this system increase with the initial errors. 
This increase is nonlinear as is to be expected since the computation of the 
crossrange deviation from the state vector at the time of entry is nonlinear. 
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The fixed-landing-site guidance system does constrain crossrange 
adjustment^ and the rms crossrange deviation is determined mainly by errors in 
the final velocity correction and trajectory estimation errors at the time of 
this correction. The crossrange deviations for fixed- landing-site guidance 
are much more sensitive to errors in the final velocity correction than those 
for fixed time of arrival. However^ the total deviation for landing-site 
guidance is so small compared to that for fixed time of arrival that this 
increased sensitivity to the velocity correction errors is not significant. 

Dowirange .- The rms entry- range deviations for the tvo systems and the 
rms entry-range adjustment for fixed- landing-site guidance are presented in 
table YE . As in the case of perigee altitude and crossrange deviations the 
entry-range deviation for the f ixed-time-of-arrival system is much less sensi- 
tive to errors in the final velocity correction. The entry- range adjustment 
for fixed- landing-site guidance results from adjusting the entry-range angle 
and azimuth angle to compensate for deviations from the reference position. 
Hence ^ the rrnis value of the adjustment increases linearly vith the square root 
of kg. The total deviation in table VI is the root sum square of the entry- 
range deviation and the do-wnrange adjustment . Under the assimiption that the 
two components are independent random variables this quantity represents the 
total rms change from the reference entry range required of the vehicle. It 
can be seen by comparison of cases la^ 2a^ and 3a vith cases Ib^ 2b ^ and 3h 
that the contribution of the entry-range deviation^ particularly the part 
caused by velocity correction errors^ is relatively unimportant. 

Landing time .- The rms landing -time deviations for the two systems and 
the rms landing-time adjustment for fixed- landing-site guidance are given in 
table VII. Like the entry-range deviation^ the landing-time deviation is due 
principally to errors in the final velocity correction and in trajectory esti- 
mation at the time of that correction. In this case^ in contrast to the other 
landing parameters^ velocity correction errors have negligible effect on the 
fixed-landing-site guidance system. In addition^ for fixed- landing- site guid- 
ance the landing-time deviations are negligible compared to the landing time 
adjustment. Since the landing-time adjustment and the entry-range adjustment 
are closely related^ the landing-time adjustment also increases approximately 
linearly with the square root of kg. 

Effect of changing observation and veloci ty correction schedule.- The 
terminal accuracy of the fixed- landing-site guidance system with the standard 
schedule of observations and velocity corrections is adequate; however^ it was 
desired to compare the total corrective velocity requirements needed when two 
systems provide approximately the same terminal accuracy. As was pointed out 
earlier^ the perigee altitude error can be reduced by delaying the final veloc- 
ity correction^ but the magnitude of the final correction increases rapidly as 
the delay is increased. It was found empirically that the perigee altitude 
error can be reduced by making additional observations during the delay time. 

A few trials with this procedure resulted in the following revised schedule: 
Starting one hour after the last observation for the standard schedule^ eight 
observations were made of the earth followed by the third velocity correction^ 
all at half-hour intervals. This revision delayed the final correction until 
about 4.6 hours before entry. 
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When kcr = 1 vas used with the new schedule^ the terminal errors of the 
fixed- landing-site system were comparable to those of the fixed-time-of- 
arrival system for case lb (e.g.^ the rms perigee altitude error was reduced 
from 2,11 to 1.57 km). This improved accuracy required an increase from 0.26 
to 0.53 m/sec in the rms value of the final velocity correction. The total 
rms corrective velocity for fixed-landing-site guidance ranged from about 
62 percent of f ixed-tirae-of-arrival guidance with the standard schedule in 
case lb to 66 percent in case 3h. Thus^ even with the increase in the final 
correction^ the fixed- landing- site guidance system still saves a substantial 
portion of the midcourse fuel. 

Computer requirements .- The fixed-time-of-arrival system used the same 
computer program as was used for the data in reference 2. A simple modifica- 
tion of this program was required to obtain data for the fixed-landing-site 
guidance system. The increase in program size was about equivalent to the 
amount required to compute the multiplying factor (l - uuT) and multiply the 
matrix (M -l) by it. Wo noticeable change in computation time resulted. 

The initial conditions required for the fixed-landing-site guidance system are 
the same as those for fixed-time-of-arrival system with the addition of the 
matrix ^1;^ discussed in appendix A. 


CONCLUDING REMARKS 


The fixed- landing-site guidance system has been shown to require only a 
little over half as much corrective velocity as the fixed-time-of-arrival sys- 
tem in the presence of injection errors considered to be realistic. This 
ratio increases to about 2/3 for an order of magnitude increase in injection 
errors^ but because the total corrective velocity also increases^ the actual 
saving in fuel is greater. 

The fixed- landing- site system is more sensitive to errors in the final 
velocity correction than the fixed-time-of-arrival system^ but this increased 
sensitivity will not signif icantly affect the objective of making the desired 
landing. It is possible to improve the terminal accuracy to that of the 
fixed-time-of-arrival system by revision of the observation and velocity cor- 
rection schedule. This improvement increases the corrective velocity required 
to about 2/3 of that of the fixed-time-of-arrival system. An alternate method 
would be to use a rocket engine for the midcourse correction having a smaller 
cutoff error than the one assumed in this study. This approach would not 
degrade the total corrective velocity requirement^ but would reduce the error 
caused by errors in making the final velocity correction. 

The errors assumed in the study are considered to be realistic^ but it 
must be noted that the relative rms corrective velocity required by the two 
systems is strongly dependent on the distribution of errors. The relative 
merits of fixed- landing- site guidance can only be determined on the basis of 
the type of observed data and the statistical distribution of injection errors^ 
observation errors^ and the velocity correction errors. 
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The computations required are of about the same complexity for the two 
systems. However^ the computation of the input data is slightly more compli- 
cated for fixed- landing-site guidance. 

The example presented here was an application to return to earth, from the 
moon^ but the basic method can be applied to direct landing on other rotating 
bodies . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 3 , I965 
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APPEKDIX A 


DERIVATION OF THE FIXED-LANDING~SITE GUIDANCE EQUATION 


The first step in obtaining the linearized fixed- landing-site guidance 
equation is the derivation of the relationships given in generalized form in 
equations (l) . For this purpose the angle from vacuum perigee to landing is 
defined as cpp. It is assumed that the angle from entry to vacuum perigee is 
fixed at 12^ (cpp = cp - 12^) and that it is traversed in 122 seconds. Thus^ 
from equation ( 2 ) the time of vacuum perigee is 


tp = tL - 0.0093cpp - 0 . 00214-0 


(Al) 


We now define three coordinate frames, illustrated in sketch (c) by three 



sets of unit vectors: (l) Up along 

the radius vector to vacuum perigee^ 
Uy along the perigee velocity vec- 
tor^ and uj 2 normal to the orbital 
plane; ( 2 ) u^ along the radius vec- 
tor to the location of the landing 
site at the time of landing^ Ux 
eastward from the landing site^ and 
Us northward from the landing site 
( an auxiliary unit vector to this 
set is U 3 which is normal to u-|^ 
and in the orbit plane); ( 3 )_an 
inertial Cartesian set with_ I 
toward the vernal equinox^ K north- 
ward along the earth ^s axis^ and J 
completing a right-handed orthogonal 
system. 


Sketch (c) _ unit vectors up^ u^, and 

U 3 lie in the orbital plane and it 

can be seen from the sketch that Up can be resolved into components cos cpp 
and sin cpp along ul anU U 3 ^ respectively. Now Ux^ Us^ and U 3 lie in a 
plane normal to u-j^ (tangent to the earth^s surface) so that the component of 
Up along U 3 may be further resolved into components along Ux and Us to 
give 


Up = up cos cpp + Ui sin cpp sin a + Us cos cpp cos a (A2) 

Similarly^ it can be shown that 

Uy = -Uj^ sin cpp + Ux cos cpp sin a + u^ cos cpp cos a (A3) 

From sketch (c) it can be seen that the Uj^, Ux^ u^ system is related to 
the Cartesian system by the transformation. 
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UL 


cos Dl cos RAl 

cos Dl sin RAl 

sin Dl 

i 



-sin RAl 

cos RAl 

0 

j 



.-sin Dl cos RAl 

-sin It sin RAl 

cos Dl. 



Equations (A2)^ (A3)^ and (a 4) can be combined to give 

Up = l(cos I>L cos RAl cos cpp - sin RAl sin 9p sin a 

- sin Dq, cos RA^ sin cpp cos a) + j(cos Dl sin RAl cos cpp 
+ cos RAl sin 9 p sin a - sin Dl sin RAl sin cpp cos a) 

+ K(sin Dl cos cpp + cos Dl sin cpp) 


(Ah) 


(A5) 


Uy = I (-COS Dl 9p “ sin RAl cos cpp sin a 

- sin Dl cos RAl cos cpp cos a) + J(-cos Dl sin RAl sin cpp 
+ cos RAl cos cpp sin a - sin Dl sin RAl cos cpp cos a) 

+ K(-sin Dl sin cpp + cos Dl cos cpp cos a) (a6) 

Rewriting equations (1) ve have 

Rp = RpUp 


which together with equations (A5) and (a6) relate perigee position and veloc- 
ity to the arrival parameters. Note that the assumption that Up and uj^ ^oth 
lie in the orbital plane implicitly restricts the incoming orbit to one that 
requires no crossrange adjustment. 

We now wish to linearize the problem by considering first-order 
perturbations about a reference trajectory subject to the constraint that the 

perturbed trajectories have the same 
values of Rp and Dp as the reference. 
Differentiating equations (AT) will 
result in a matrix of the first partial 
derivatives of perigee position and 
velocity with respect to the arrival 
parameters . However, we wish to obtain 
first partial derivatives of position 
and velocity at the time of reference 
perigee. The distinction between these 
quantities is illustrated with the aid 
of sketch (d). 
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Here Rp and Vp are position and velocity vectors of perigee on the 
perti^bed trajectory which occurs at some time ta* On the other hand^ R(tp) 
and V(tp) are the position and velocity vectors on the perturbed trajectory at 
the time^ tp^ when reference perigee occurs. If ta is nearly equal to tp^ 
then we can write 

R(t ) = Rp - V (ta - tp) = Rp - Vp(6ti - 0.00933&p) 

(A8) 

V(tp) = Vp - Vp(ta - tp) = Vp - Vp(6ti - 0.00933&P) 

Define \|/p as the matrix of first partial derivatives of position and 
velocity at the time of reference perigee with respect to the arrival param- 
eters. We now compute ^|/p by differentiating equations (a 8)^ but we must 
recognize that the terms involving (tg^ “ "*^p) already in the form of first- 

order perturbations. Therefore^ 



The computation of the partial derivatives in equation (A9) is simplified 
by expressing Rp and Vp as in equations (l) . 

Rp = RpUp 

Vp = VpU^ 



IT 



A change in the angle a 
Ul so that 


represents a rotation of the orbital plane about 



S<x 


= Rp X Ul 


^ = V — 



X Ul 




> 


J 


(A12) 


A change in the time of landing rotates the orbital plane about the earth’s 
axis. If cjo is the angular velocity of the earth then 


BtL ^ ^■'^L 


^ 

Str. ^ 


Rp(wK X Up) 


Vp(u)K X u^) 


(A 13 ) 


Finally^ 


avj 




0 


Uv 




(AlU) 


Substituting equation's (All) through (A14) in equation (A3) with the 
unit vectors in terms of Rp and Vp gives 


^n = 


Ar 


V^P 


£ + 0.0093) V Rp X Uy (K X Rp)-Vp 




V p 

+ 0.0093 V. 


Vp X 


(K X Vp)-Vp 






This expression can be verified by straightforward differentiation of 
equations (A5) and (a 6) followed by considerable algebraic manipulation. 

We will now use to derive the linearized fixed- landing-site guidance 

equation. For this purpose we define r and v as small vector deviations 
from the reference position and velocity and x as the 6x1 state vector 
given by 


X = 


r 


(A15) 
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The state vectors at different times are related by the state transition 
matrix 0 as follows : 


(Al6) 


x(ta) = ^>(t2^ti)x(ti) 

If ts is the present time_, t^ and ti the time of reference perigee^ 
tp^ then equation (A16) becomes 

x(t) = 0(t,tp)x(tp) (ait) 

Let ^ be defined as a 4 x 1 toatrix of deviations from the reference 
arrival parameters . Then 



" V 


substituting for 

x(tp) in equation (AI 7 ) 



x(t) = «>(t,tp)\|fp^ 

(AI 8 ) 

or 

X = 

(AI 9 ) 

where 

\|; = 4>(t,tp)l/p 

(A20) 


Equation (AI 9 ) may be -written in partitioned form as 

^2 

u 


r' 


^1 Pi 

V* 


^2 Vz 


(A21) 


Here and ^ 2 . 3x3 matrices, Pi and "ps ar® 3x1 matrices and the 

indicate deviations from the reference values of the arrival parameters . The 
primes have been added to r and v to indicate that they are the deviations 
from the reference position and velocity -which would be produced by changes 
of ^ 1 , 1 ^ 2 } ^ 3 } ^4 ■fcli® arrival parameters. 


Expansion of equation (A21) gives 

r: 


= 


Ud 


+ Pl^4 


(A22) 
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(A23) 


If is taken to be the actual deviation of the vehicle's position from the 
reference trajectoiy^ the changes in the first three arrival parameters neces- 
sary to produce a trajectory from the present position to the landing site are 


^2 = r - 


(A24) 


The velocity deviation of this trajectory from that of the reference trajec- 
tory is found by substituting equation (A2U) into (A23) 

v» = Pi - Pa)U 

The total difference between the desired velocity and the present velocity is 

v„ = v» - V 


or if 


M = \|f2't'j 


_ Pi ~ P2) “(Mpi “ P2) 


^2^1 ^ Pi “ P2 


then 


V = Mr + ku^4 - V 


(A25) 


It can be seen from equation (A25) that v^ is not uniquely determined 
but is a function of ^4. It is desired to choose ^4 so as to minimize 
I Vq I ^ and this optimum is V(f[ in equation (3) 

+ Ku^ 4 - v)^(Mr + ku ^4 - v) 

^ = 2k^^4 + 2ku^(Mr) - ku% 
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(A26) 


If this derivative is set equal to zero, the optimum 


u 


-(u^iylr) “ u^ 


^4 is found to be 


Substituting in equation (A5) gives 

— T -- _T— — 

= jyir - uu Nr + uu V - V 


or 
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APPENDIX B 


ERROR MODEL FOR THE DIGITAL COMPUTER SIMULATION 


The error model assiimed for mechanizing the velocity correction can be 
explained briefly 'with the aid of sketch (e) in which v^^^ is the desired 

velocity correction while Vg^ is 
the correction actually made. 

The directional displacement 
between v^ and Vg^ can be speci- 
fied by the angles e and 6. The 
angle^ between the two vectors 
may be either positive or nega- 
tive and has Gaussian distribu- 
tion^ zero mean^ and 1^ standard 
deviation. The angle^ gives 
the orientation of the plane of 
the two vectors with respect to 
some fixed reference and has uniform distribution between zero and at. There 
is also an error in jvg| which has two components. There is a cutoff error 
having Gaussian distribution^ zero mean ^ and a standard deviation of 0.1 m/sec. 
The other component due to thrust level also has Gaussian distribution with 
zero mean^ but its standard deviation is 1 percent of bdl- 



The relative significance of the different components of the velocity 
correction error can be better explained with the aid of sketch (f) showing 



Sketch (f) 

Likewise^ Vi is approximately equal to 


Vi + 


the plane of v^^^ and Va* Th^ error ^ 
Vg^ in the velocity correction can be 
resolved into two components^ v^^ in 
the direction of V(i and v^^ normal to 
Since e is very small a good 
approximation to v^ is given by 

V2 VdG 

the error ^ y^^ in Vg so that 
(O.Ol Y^)^ 


Since the rms value of e is 1^ (about l/60 radian)^ it is seen that for v^^^ 
of order 1 m/sec^ v^ is nearly an order of magnitude smaller than v^^ and 
the rms value of Vi consists almost entirely of the cutoff error. On the 
other hand^ if v^^^ is of order 10 m/sec^ the irnis value of v^ is about 
0.17 m/sec while that of Vi is about 0.1^ m/sec. 


The method of calculating the covariance matrix^ S^ of velocity correction 
mechanization errors from the error model just described is given in 
reference 5« 
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The errors in measuring the individual Cartesian components of velocity 
correction were assumed to be random and independent with Gaussian distribution^ 
zero mean^ and standard deviation of these measirrements is 

Sm = 

As mentioned in the text it was assimed that certain celestial 
observations (angles) are made for the purpose of estimating trajectory. These 
measured angles consist of the subtense angle of either earth or moon together 
with the azimuth and elevation angles of center of the same body with respect 
to some inertial reference. 

The errors in the observations were assumed to be random with Gaussian 
distribution and zero mean. Also^ the errors in different angles and measure- 
ments made at different times were assumed to be uncorrelated. The standard 
deviation was the same for the three angles measured and is assumed to be given 
in seconds of arc by 

a = 

where 7 is half the subtended angle of the body observed. These statistical 
characteristics result in the covariance matrix^ of observation errors being 
a diagonal 3x3 matrix 

Q = a^I 
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APPMDIX G 


COMPUTATION OF STATISTICAL INFORMATION FOR ORBITAL PHASE OF FLIGHT 


During the orbital phase of the flight the statistical performance of both 
guidance systems was evaluated by means of linear theory. This evaluation 
required the knowledge of three covariance matrices: 

(1) The covariance matrix^ of position and velocity deviations from 
the reference trajectory^ 

(2) The covariance matrix^ P^ of errors in estimating those quantities^ 

(3) The covariance matrix^ C^ of velocity correction actually made. 

This third matrix is separated into two parts so that the penalties due to 
velocity correction errors may be evaluated. A supplementary covariance matrix 
is needed^ for reasons to be discussed later^ to account for adjustments made 
in entry range and landing time by the landing-site guidance system. 

The covariance matrices and Q of velocity correction mechaniza- 

tion errors^ velocity correction measurement errors^ and errors in celestial 
observations are used in calculating P^ and C. The sources of these 
matrices are given in appendix B. 

Initially^ P and W are both set equal to the covariance matrix of 
injection errors^ and they must be updated as the flight proceeds. This 
updating which takes place either as a result of the passage of time or because 
of operations (observations^ velocity corrections^ or termination of the flight) 
performed by the system proceeds as follows: 

Let ^ be the transition matrix giving deviations from the reference 
position and velocity at time t due to deviations at some earlier time. Then 
between operations^ P and W are updated by 

P = 

¥ = $¥o<I>^ 

If H is the matrix of partial derivatives of the observed angles with 
respect to the Cartesian components of position and velocity^ then after an 
observation 

P = Pq - PoH^LhPoHT + Q]"Vo (C 2 ) 

while W remains unchanged. The derivation of this relationship is quite 
complex and will not he presented here. (See ref. k, eqs . (13) and (l 4 ).) 
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The only estimation errors incurred in making a velocity correction are 
those in measuring the correction actually made. Hence ^ after the correction 


P = Pq + 


(C3) 


It is shown in reference 2 that straightforward computation of E(xxT) 
after a velocity correction gives 


W = 


I 0 

(G 2 + l) 


(Wo - Pq) 


I 0 

Gi (Ga + I) 


+ Po + S 


(C4) 


Where I is a 3 X 3 unit matrix^ 0 represents a 3 x 3 null matrix^ and 

G = (Gi G 2 ) 


The significance of the different terras on the right side of equation (C4) can 
be explained by the use of equations (Cl) to update W to the time of refer- 
ence perigee . Then 


W(tp) = Wi + Wa 


where by definition 


W-L = ^ 


Gi (G2 + l) 


(Wo - Po) 




0 


Gl (G2 + l) I 


and 


Wa = 0>(Po + 


It is also shown in reference 2 that (W-P) is the covariance matrix of 
the estimated deviations from reference position and velocity. Since the 
velocity correction would cause the guidance requirements to be satisfied if 
the estimate were correct^ Wi represents a distribution of deviations from 
the reference trajectory for which the guidance requirements are still satis- 
fied. On the other hand W 2 represents deviations^ from members of this fam- 
ily of satisfactory trajectories^ which are caused by errors in trajectory 
estimation at the time of the velocity correction and by errors in making the 
correction. 


From this analysis it appears that the first term on the right side 
of equation (c4) could be eliminated^ with no loss of significant information^ 
to give 


w = Pq + S (C5) 

However, the information contained in the term eliminated from equation (c4) is 
necessary for computing the statistics of the landing parameters. The above 
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reasoning can be confirmed analytically for the case of fixed-time-of-arrival 
guidance as follows: 

For the fixed-time-of-arrival guidance equation stated in equation (2) 


G = —0 2 ^ 1 

G2 = -I 


(c6) 


Here I is the unit matrix and and ^2 arise from partitioning as 
follows : 


0) = 


I ^ 1 ^ 2 . 

^3 ^4 

Substitution from equations (c6) and (C7) shows that 


0 


I 0 

Gi (^2 + 

It is shown in reference 2 that 


d) 


-1 


* 4 ’’ 

- 43 ^ 


and this relationship can he used to show that 


0 


(O 3 ” d>4d>2 ^ 1 ) 0 j 


(CT) 


(C8) 


(C9) 


($3 — ^ j_) — (“"fH 


(CIO) 


Therefore^ for the fixed-time-of- arrival guidance equation 



0 


0 d >2 

Wi = 


(Wo - Po) 



($2"^)^ 0 

0 

1 


(Cll) 


Hence^ only the lower right 3x3^ the velocity components, of are not 

zero. These are the deviations in velocity at the time of reference perigee 
which must be accepted in order to have zero position deviation. Since the 
landing parameters are functions of perigee velocity as well as position, it 
is necessary to retain the information contained in This information was 

retained in the case of f ixed-time-of-arrival system by using equations (Cl) 
to update W. (Since, to a first-order approximation, d >2 is proportional to 
the time over which it is computed, the effects of Wi on the landing 
parameters will be small.) 


In the case of fixed-landing-site guidance it would be very difficult to 
reduce mathematically to a simple form. However, as in the case of 
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f ixed-time-of-arrival system^ will make no contribution to the deviations 
which the guidance law requires to be zero. For this system these are devia- 
tions in perigee altitude^ declination of the landing site_, and crossrange 
adjustment. In this case the information lost by the elimination of will 

be the statistical distribution of the adjustments in entry-range angle ^ land- 
ing time^ and perigee velocity made by the guidance system in order to achieve 
a satisfactory trajectory. It was desired to know the first two of these three 
quantities which are defined in the text as downrange adjustment and landing- 
time adjustment^ respectively^ but preliminary results using equations (Cl) 
with the fixed-landing-site guidance system showed that round-off errors cause 
unacceptable inaccuracy for large initial values of W, For this reason equa- 
tion (C5) was used with this system and the statistics of the downrange adjust- 
ment and landing-time adjustment were computed separately as described in the 
following paragraphs. 

The generalized inverse of the matrix ij; in equation (AI 9 ) (derived in 
appendix E) can be used to compute these statistics. This inverse will be 
written so that 

^ = il/*x (C12) 

where ^ is the result of applying the guidance equation. 

At the time of any given velocity correction 

Note that this is the covariance matrix of ^ due to a specific velocity cor- 
rection^ and W is evaluated immediately before that correction. Under the 
assumption that the various velocity corrections are uncorrelated^ the covari- 
ance matrix for the total change in arrival parameters can be found by summing 
the individual covariance matrices. 

The mean-square velocity correction for the nth correction time is 
given by the trace of the covariance matrix. 

Cn = G(Wn - Pn)G^ + (Clif) 

In equation (Cl4) and are evaluated just before the correction. If 
is obtained by updating W from the time of the previous correction then 

Cn = G($Wn-i$^ - Pn)GT + 

Substitution for 'Wn_j_ from equation (C4) yields 

Cn = G4)(Wn_^ - Pn-i)5>^G^ + G<|)Pn_i3>^G^ + CM'Sn-ii’^G^ + Sn (Cl6) 
where and are evaluated immediately before the n-lst correction. 


(C15) 


(C13) 


2T 



(Also note that the first term on the right side of equation (Cl6) should he 
zero.) The trace of the matrix gives the mean-square value of 

that portion of the nth velocity correction resulting from errors in making 
the previous correction. 
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APPENDIX D 


COMPUTATION OF STATISTICS OF THE LANDING CONDITIONS 


The application of the Monte Carlo method to the computation of the 
landing conditions is as follows: 

1. Starting immediately after the final velocity correction^ two-body 
formulas were used to compute the perigee radius of the reference 
trajectory and its position and velocity vectors at atmospheric entry. 

2. The formulas derived below were used to compute the reference landing 
conditions . 

3 . A set of random position and velocity deviations from the reference 
trajectory was obtained from the statistical distribution described 
by W (ref. 6). These deviations were added to the reference posi- 
tion and velocity to obtain a perturbed position and velocity. 

4. Steps 1 and 2 were repeated for the perturbed trajectory and the 
differences between perturbed and reference landing conditions were 
computed. 

5. A sample of 1000 different perturbed trajectories was found to give 
sufficiently accurate values of the statistical quantities of 
interest . 

Standard two-body formulas were used for Step 1 and will not be presented 
here. The entry position^ velocity^ and time were then used to compute the 
landing conditions and the formulas used for this computation are derived in 
the following paragraphs . 

The relationship of the unit vector u^^ along the radius vector to the 
landing site to the Cartesian coordinate system is shown in appendix A to be 

UL = I cos Dl cos RAl + J cos Dx, sin RAj, + K sin 

Similarly^ the unit vector u;e along the radius vector to the entry point is 
given by 


U;g = I cos D;g cos RA;g + J COS Dg sin RA;g + K sin Dg 


Since 


cos cp = U;g • ul 


it follows that 

cos cp = cos Dl cos Dg cos (RAl - RAg) -f* sin Dl sin Dg 


(D2) 


(D3) 
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If RA-le is the right ascension of the landing site at the time of entry, 
then 

RAl = RAle + ‘^tp (d^) 

where w is the earth’s angular velocity and tp> is the time of flight from 
entry to landing. Substitution in equation (D3) gives 

cos cp = cos Dg cos Dl cos (RAt.-r; + wtp - RAe) (D5) 

and from equation (2) 


tp = 0.00933CP + 0.00254 


Equations (2) and (D5) must be solved by iteration (see ref. l) to find cp and 
tp. The time of landing, t^, is 


K 



Sketch (g) 


The computation of the crossrange 
deviation can be explained with the aid of 
sketch (g). 

If the vehicle were allowed to continue 
in its incoming orbital plane the ^'actual 
orbital plane" in the sketchy it would cross 
the landing site track at point Li. At the 
same time the landing site would be at point 
Lq* If the vehicle were on a satisfactory 
trajectory^ it would pass through the entry 
point in the "desired orbital plane" and 
arrive at L2 at the same time as the land- 
ing site. The crossrange deviation is 
defined as the normal arc^ CR^ from the 
landing site to the actual orbital plane. 


Let ng^ be a unit vector normal to the actual orbital plane and n^ a 
unit vector normal to the desired plane. Then 


_ Re X Vg 

X VeI 

and 

^ 

^d - -rz ITT 

The relative inclination angle, i, can be obtained from 
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and from spherical trigonometry 


cos i = 


tan (CR) = tan i sin cp 

Multiplying by a constant factor converts CR from angular units to 
kilometers . 
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APPENDIX E 


DERIVATION OF THE GENERALIZED INVERSE OF i|r 


It is shown in reference 7 that every rectangular matrix. A, has a 
generalized inverse. A*, which is unique. Furthermore, it is a property of 
A* that if 


b = Ay (El) 

has a solution, then that solution is given by 

y = A*b (E2) 


The guidance law was derived to insure that for every state vector x 
there exists a unique Therefore, it is known that 


5 


(E3) 


This fact will be used to solve for \lr*. 

Equation (E3) can be written in partitioned form as 


is 

ia 

u 

where ijfi* and ts* 3x3 


) 

1 

* 

r 

/- *.,T 
(Pi ) 

(ps*f 

V 


matrices. From equation (A26) 


(e4) 



SO that 


k 

k 


(E5) 


If the optimum ^4 from equation (A26) is substituted into equation (A22) 
(the prime is omitted from r) then 


r = 


^3 


- — pi(u^Mr - u^v) 
k 
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so that 


1 1 

UsJ 


From equation (e6) 


= ^(I + 7 Piu%)r - \)/i ^ -^lU^ 

k k 


^1* = ^ PiU^M) 

ta* = - Piu^ 

k 


(e6) 


Substitution from equations (ET) and (E5) into equation (e 4) gives 


(ET) 


\1/* = 


i piu%) - i 


1 ttT 


— (u%) — u 

k k 


It can be shown that 


= I 


where I is a ^ X 'unit matrix and 

= I 

where I is a 6 x 6 mit matrix . 
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TABLE I.- COVARIANCE MATRIX OF INJECTION ERRORS FOR STANDARD DISTRIBUTION 



X, 

km 

1 

X, ! 

km 1 

z, 

km 

X, 

km/sec 

km/sec 

• 

km/sec 

X 

0.50443614 

-0.1053457 

-0.16956842 

-0.73958537x10"^ 

-0.47001117x10"® 

-0.22777452x10"® 

Y 

-.10534570 

. 36457942 x 10 "^ 

.30887450x10"^ 

. 17144754 x 10 "® 

.15240230x10"® 

.59077529x10"^ 

Z 

-.16956842 

.30887450x10"^ 

.84712863x10"^ 

.29243600x10"^ 

.17300981x10"® 

.13020790x10"^ 

• 

X 

-.73958537x10“^ 

.17144755x10"^ 

.29243600x10"^ 

.22371290x10"® 

-.11326999x10"® 

.25554670x10"® 

• 

Y 

-.47001117x10"® 

.15240228x10"® 

.17300981x10"® 

-.11326999x10"® 

.18436915x10"® 

.25527313x10"® 

• 

Z 

.59077530x10"^ 

.59077530x10“^ 

.13020790x10"® 

.25554670x10"® 

.25526313x10"® 

.15039^92x10"^ 


vn 



OJ 

ON 


TABLE II.- BMB 'VELOCITI CORRECTIONS 



Fixed-time-of-arrival 
correction no. 

Fixed-landing-site 
correction no. 

Ratio of 
total for 
F * L • S • to 
that for 

F.T.A. 

Case 

ka 

Velocity 

correction 

errors 

Injection 

errors 

m/ sec 

2, 

m/sec 

3 , 

m/ sec 

Totals 
in/ sec 

1. 

m/ sec 

2, 

m/ sec 

3 . 

m/sec 

Total ^ 
m/ sec 

la 

1 

None 

standard 

1.37 

1.08 

0.72 

3.17 

0.94 

0.55 

0 .l 4 

1.63 

0.52 

2 a 

10 

None 

Standard 

4.58 

1.32 

.75 

6.65 

3.15 

.66 

.15 

3.96 

.59 

3 a 

100 

None 

Standard 

l 4 . 6 l 

1.64 

.84 

17.09 

10.04 

.77 

.15 

10.96 

.64 

lb 

1 

Standard 

Standard 

1.37 

1.11 

.83 

3.31 

. .94 

.58 

.26 

1.78 

.54 

2 b 

10 

Standard 

Standard 

■ 4.58 

1.36 

.86 

6.80 

3.15 

.69 

.27 

4.11 

.60 

3 b 

100 

Standard 

Standard 

l 4 . 6 l 

1.76 

.94 

' 17.31 

10.04 

.83 

.27 

11.14 

.64 

1 

1 

Standard 

Spherical 

2.79 

1.47 

.93 

5.19 

2.68 

.71 

.27 

3 .66 

.70 


100 

Standard 

Spherical 

28.62 

2.47 

.99 

32.08 

27.22 

1.15 

.27 

28.64 

.89 


TABLE III.- CORRECTIONS DUE TO ERRORS IN PREVIOUS CORRECTION 



Fixed time of arrival 

Fixed landing site 

Case 

Second correction, 
m/sec 

Third correction, 
m/sec 

Second correction, 
m/sec 

Third correction 
m/sec 

lb 

0.24 

0.40 

0.16 

0.20 

2b 

.31 

. 4 o 

.18 

.20 

3 b 

.63 

.41 

.30 

.20 



TABLE IV.- RMS PERIGEE ALTITUDE DEVIATIONS 


Case 

la 

2a 

3a 

lb 

2b 

3b 


Rms perigee altitude deviations^ 
'km 


Fixed time of arrival 

1.51 

1.54 

1.50 

1.56 

1.58 


Fixed landing site 

1.49 

1.51 

2.11 

2.17 

2.20 


TABLE V.- RMS CROSSRAWGE DEVIATION 



Crossrange deviation^ 
km 

Case 

Fixed time of arrival 

Fixed landing site 

la 

3.82 

0.59 

2a 

5.61 

.61 

3a 

12.26 

.62 

lb 

3.85 

.92 

2b 

5.67 

.93 

3b 

12.37 

.93 
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aero7iautical and space activities of the United States shall be 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.*^ 

— National Aeronautics and Space Act of 1958 
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